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A numerical simulation is presented here of the evolution of initial data of the kind that was 
conjectured by Hertog, Horowitz and Maeda [J to be a violation of cosmic censorship. That initial 
data is essentially a thick domain wall connecting two regions of anti-deSitter space. The initial 
data has a free parameter that is the initial size of the wall. The simulation shows no violation of 
cosmic censorship, but rather the formation of a small black hole. The simulation described here 
is for a moderate wall size and leaves open the possibility that cosmic censorship might be violated 
for larger walls. 
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j^jy One of the outstanding questions in general relativity is that of cosmic censorship: the question of whether the 
singularities formed in gravitational collapse are hidden inside black holes. One thing that makes this issue difficult 
is that it is not clear what we might expect to be true, ft is fairly straightforward to find spacetimes with naked 
singularities; however, it can certainly be argued that those spacetimes are not generic (as in the case of critical 
collapse 2]) or that the matter is not "physically reasonable" (as in the case of dust). The question is what is to 
be meant by "physically reasonable" and whether cosmic censorship holds if the condition of physically reasonable 
matter (as well as genericity) is imposed. One might think that pointwise energy conditions would be a condition that 
would be satisfied by physically reasonable matter. However, Hertog, Horowitz and Maedap] point out that matter 
can violate pointwise energy conditions and still give rise to a positive mass theorem. They argue that such matter 
■ should be considered physically reasonable; but that it gives rise to naked singularities. In reference Q] initial data 
is generated and it is argued that the evolution of that data gives rise to a singularity that cannot be hidden inside 
i a black hole. Such a singularity would have to either (i) be visible to an observer at infinity or (ii) itself extend to 
infinity. In reference^] it is argued that possibility (i) is more likely; however, a recent theorem of DafermosQ rules 
out possibility (i). Thus if the arguments of reference^] are correct then their initial data evolves to a singularity 
that extends out to infinity. The system is gravity coupled to a scalar field with a potential that has two minima 
below zero. The initial data is essentially a thick domain wall that interpolates between the true vacuum on the inside 
and the false vacuum on the outside. The arguments of reference yj are somewhat heuristic, so it is worthwhile to 
evolve their initial data to see whether a naked singularity does indeed form. Alcubierre et al perform numerical 
simulations to evolve an initial data set for a system that has some similarities to the system considered in reference[l|. 
Here the system is also gravity coupled to a scalar field with a potential with two minima; and the initial data is a 
thick domain wall that interpolates between the two vacuua. However, in this case one minimum of the potential is 
at zero and it seems unlikely that the positive mass theorem holds for this system. The results of are that the 
domain wall accelerates outward, approaching the speed of light. 

Since the results to be obtained may depend crucially on the type of matter used, it is important to simulate the 
system described in reference^ rather than a different system that may or may not be analogous. In this work we 
perform a preliminary version of such simulations. The work is preliminary for the following reason: the argument of 
reference[l| says that the singularity cannot be hidden inside a horizon if the initial wall is sufficiently large. Here, 
sufficiently large means a radius greater than about 600 in units of the anti-deSitter radius of curvature. Simulating 
such a large wall is numerically very challenging. Instead, this simulation is for a wall of the more moderate initial size 
of about 7. The system and data, as well as the numerical methods arc described in section II. Results are presented 
in section III and conclusions in section IV. 
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II. EQUATIONS AND NUMERICAL METHODS 



The system to be studied is a spherically symmetric scalar field <f> with a potential V. The appropriate equations 
are therefore the Einstein-scalar equations: 

G ab = V Q 0V 6 0- gab{\V c $V c c^ + V) (1) 
dV 

V B V* = ^ (2) 

(Here we are using units where 8irG =1). We use polar-radial coordinates for the metric which puts it in the form 

ds 2 = ~a 2 dt 2 + a 2 dr 2 + r 2 (d6 2 + sm 2 6d(p 2 ) (3) 



It is helpful to define the quantities X and Y given by 



Y = - f (5) 



a 
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from which it follows that equation of motion for <f> is 



°Y (6) 



dt a 

From equation J5J it follows that the equation of motion for Y is 

dY 1 d ( 2 a \ dV 

~dt = ^{ r a X )- aa dj (7) 
From equation JQl we obtain the following constraint equations for a and a 

da a(l — a 2 ) 



dr 2r 



(X 2 + Y 2 + 2a 2 V) (8) 



Aln(aa) = ^ (X 2 + F 2 ) (9) 

We now turn to numerical methods. The equations simulated are equations We use unequally spaced values 

of r. At grid point i the quantity X is calculated as Xi = {4>i+\ — 4'i~i)/( r i+i — r %-i)- For any quantity / define 
quantities /+ and /_ by /+ = (fi+i + /»)/2 and /_ = (fi + /j_i)/2. Also define the quantities X p and X m by 
X p = (4>i+i — (j)i)/(r i+ i — r%) and X m — (fa — fa-i)/(ri — The first term on the right hand side of equation Q 

is differenced as 



r z or \ a / ri — r_ \ ct+ ct_ 



X ) - n ( ^+ — - r - — *«» ) (10) 



Equations Ij6l7|l are integrated using a three step iterated Crank-Nicholson method0 with Kreiss-Oliger dissipation Q 
used in equation Q . Each iteration of the Crank-Nicholson method involves an integration of equations 1-SI!)[) which 
is done using a second order predictor-corrector method. 

At the origin, (gridpoint i = 1) smoothness requires that <fi and Y have vanishing derivative with respect to r. 
We impose this as fa — (Afa ~ <fe)/3 and correspondingly for Y, We also impose the condition a\ = 1 which is 
required by smoothness of the metric. The value of a at the origin can be freely specified, and we choose a.\ = 1. At 
large distances, the spacetime should be in the anti-deSitter false vacuum corresponding to <f> = </> va c- We apply this 
condition by imposing fay — fa ac and Y/v = where N is the last gridpoint. 

Stability of the simulations requires that the time step satisfy a Courant condition. We choose dt to be the minimum 
over all gridpoints i of ^(rj+i — rj)ai/a!i. This Courant condition is one reason for the use of unequally spaced r,. 
The presence of vacuum energy leads to a large variation in a/ a which would then lead to an extremely small time 
step if we had a uniform spacing for r. 

We now turn to a consideration of the initial data. The potential used in reference is 

V(fa) = -3 + 5O0 2 - 81</- 3 + kfa (11) 
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where k is a constant whose value will be specified later. The false vacuum is at <fi — while the true vacuum is at 
some other value <fi — 4> vac . The approach of reference to choosing initial data is the following: consider initially 
static field configurations that are in the false vacuum at r = and in the true vacuum for r > R vac for some constant 
i? vac ■ Then the contribution of the potential to the total mass is 

my = [\S> W2/2 VyHy (12) 

Jo 

Here y = r/R vac and 4>' = d<f>/dy. Define py to be the minimum over all field configurations of 2myR~ 3 c and choose as 
initial data the field configuration for which the minimum is attained. The positive energy theorem will hold provided 
that pv > ^(0vac)/3- For the present potential, this inequality is just barely satisfied for k — 45.928. Define x = Iny, 
Then minimizing py corresponds to the differential equation 

d 3 (j) _ d 2 (j) 
dx 3 dx 2 

We find the solution to this equation using a shooting method. Smoothness and the condition that the field be in 
the false vacuum at r = yields the condition that for large negative x we have <fi = ce@ x where c is a constant 
and /3 = (V 409 — 3)/2 « 8.6. We integrate equation Ijl3(l from large negative x to x — using the fourth order 
Runge-Kutta method. The value of c is found using a binary search to be the one that yields <f> = 4> vac at x = 0. 

The spacing in r is chosen as follows: choose a constant i? ma x and a coordinate r equally spaced from zero to i? max - 
Then define r by 

r = A 
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III. RESULTS 

The simulation was done in double precision on a SunBlade 2000. The number of gridpoints was 8001 and the 
values of the parameters were as follows: i? V ac = 7 and -Rmax = 16 while k was chosen so that </> vac = 0.725 which 
yields the critical value of k. The coordinate system used (equation Q) cannot evolve past the formation of a trapped 
surface, since it assumes that r is spacelike, whereas r becomes null on a trapped surface. As a trapped surface forms, 
the metric quantity a grows without bound. The simulation is run until a grows large enough to signal that a trapped 
surface is forming. Figure ^ shows the scalar field (f> at the initial time (solid line) and at a time shortly before the 
formation of a trapped surface (dotted line). This time is close to 7r/2. What happens in the evolution of <p is that 
the initial domain wall moves inward somewhat; but in addition, a wave pulse comes off the wall, moves to the center 
and forms a black hole there. Since the most interesting features of the final configuration are on a small spatial scale, 
subsequent figures will plot those features in the range < r < 0.25. Figures T2I4I show respectively </>, Y and a at the 
time shortly before the formation of a trapped surface. 

Note from figure [3] that the time derivative of the scalar field is getting quite large at the final time and at points 
close to the center. Thus a region of high curvature is forming near the center. Figure 01 indicates the formation of 
a trapped surface at r ~ 0.05. What is important here is not simply the formation of a trapped surface, but also 
that the size of the trapped surface is of the same order as the size of the region of high curvature rather than much 
smaller than the region of high curvature. This makes it likely that what is happening is the formation of a black hole 
that will contain the singularity. Note that the size of the black hole is considerably smaller than the spatial scale set 
by the initial size of the wall. 

IV. CONCLUSIONS 

The results of this simulation of a moderate size wall are quite different from what is conjectured in reference[l| to 
be the behavior of a large wall. Instead of a naked singularity one gets a small black hole. Instead of the collapse of 
a rather large homegeneous region, one has very localized and decidely non-homogeneous behavior near the center of 
spherical symmetry. The reason for this behavior follows from properties of the initial data: note that near the center 
the scalar field is proportional to r* 9 where (3 « 8.6. This large power means that the scalar field is essentially zero in 
a spherical region around the center. In contrast, the analysis of reference^] assumes that in the central region the 
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(13) 



l.lf + 0.9 In 



cosh(f - R vac ) 
cosh i? vac 



(14) 
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FIG. 1: the scalar field <f> at the initial time (solid line) and the final time (dotted line) 
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FIG. 2: the scalar field <j> at the final time 

scalar field can be modeled as homogeneous but time dependent. Since the initial scalar field in the central region is 
negligible, subsequent behavior of the scalar field in the central region must come from modes that propagate inwards, 
in this case from the domain wall. As these modes propagate inwards, their amplitude grows and their wavelength 
shortens until finally they form a small black hole at the center. 

Thus the results of this simulation do not provide a counterexample to cosmic censorship. Note, however that this 
simulation is done with a moderate value of i? V ac (-Rvac = 7) while the argument of reference Q claims that violations 
of cosmic censorship occur for sufficiently large R vac (R vac > 600). Thus it is still possible that this model violates 
cosmic censorship but not for moderate values of i? V ac like the one used in this simulation. Put another way, the 
arguments of reference^] do not lead one to expect the formation of a naked singularity from the i? V ac = 7 form of 
their initial data. Therefore the fact that this simulation finds black hole formation does not contradict the arguments 
of reference lj. Nonetheless, the details of the collapse process tend to cast doubt on some parts of the argument of 
reference[l|. In particular, regardless of the size of the wall, the initial data for the scalar field behaves like r 13 near 
the center for (3 w 8.6. Thus the initial data in the central region does not look like anti-de Sitter space with a small 
homogeneous perturbation. However, the argument of reference^ depends on treating the evolution of the central 
region as the evolution of a homogeneous space. Thus there is reason to doubt this part of the argument. To resolve 
this issue one would need to do a simulation with large i? V ac- Since in anti-de Sitter space a and a -1 grow linearly 
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FIG. 3: Y at the final time 
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FIG. 4: a at the final time 

with r, large values of R vac involve large values of a and a~ 1 . It is therefore likely that the numerical method used 
in this work is not sufficiently robust to treat large values of i? vac and that some other numerical method will have to 
be used. This problem is under study. 

Finally, note that the authors of reference have proposed another model0 that they claim leads to violations of 
cosmic censorship. However, numerical simulations of that model done by Gutperle and Kraus |£j indicate that that 
model leads to the formation of black holes rather than naked singularities. 
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